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SUMMARY 

The  orbit  of  the  satellite  1970-97B,  the  rocket  of  Cosmos  378,  at 
inclination  74°,  has  been  determined  at  18  epochs  between  April  and  August 
1971,  when  the  effects  of  14th-order  resonance  with  the  Earth's  gravita¬ 
tional  field  were  appreciable.  The  orbits  were  determined  with  the  PROP  6 
program  from  Hewitt  camera,  kinetheodolite,  US  Navy  and  visual  observa¬ 
tions,  and  an  average  accuracy  of  90  m  in  perigee  distance  was  achieved, 
despite  the  low  perigee  height  (near  230  km)  and  the  consequent  high  drag. 

The  orbits,  together  with  13  previously  evaluated,  have  been  analysed 
to  reveal  the  effects  of  the  14th-order  resonance  and  to  evaluate  six 
lumped  geopotential  harmonics  of  order  14.  Because  the  orbit  passed 
through  resonance  rapidly,  the  values  are  not  as  accurate  as  those  from 
slow  resonances;  but  they  are  more  accurate  than  any  others  available  for 
an  inclination  near  74 #,  and  have  proved  their  worth  in  a  recent  determina¬ 
tion  of  individual  14th-order  coefficients. 
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INTRODUCTION 


Cosmos  378  rocket1,  designated  1970-97B,  was  launched  on  1970  November  17  into  an 
orbit  inclined  at  74.0°  to  the  equator,  with  perigee  and  apogee  heights  of  233  and 
1730  km  respectively,  and  an  orbital  period  of  105  min.  It  decayed  on  1972  September  30 
after  633  days  in  orbit. 

.  2 

The  orbit  has  previously  been  determined  between  1971  February  and  decay  by  Hiller  . 

The  satellite  experienced  !4th-order  resonance  in  1971  June,  and  Hiller  analysed  the  sub¬ 
sequent  orbits  to  obtain  values  for  atmospheric  rotation  rate  and  density  scale  height. 

In  this  Report  we  analyse  the  inclination  and  eccentricity  of  the  orbits  that  are  affected 
by  14th-order  resonance. 

In  evaluating  individual  geopotential  coefficients  of  order  14  from  satellite 
resonances,  it  is  essential  to  have  values  of  lumped  harmonics  from  satellites  at  a  wide 
variety  of  inclinations,  and  previously  there  have  been  no  accurate  values  of  lumped 
harmonics  available  for  inclinations  between  70°  and  81°.  So  it  was  decided  that  analysis 
of  I970-97B  (inclination  74°)  would  be  worth  trying,  even  though  it  was  obvious  from  the 
start  that  the  values  obtained  would  not  be  particularly  accurate,  because  of  the  very 
rapid  decay  rate  resulting  from  the  low  perigee. 

A  total  of  18  new  orbits  has  been  determined  from  observations  with  the  PROP  6  pro- 
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gram  ,  for  dates  between  1971  April  3  and  August  17.  These  replace  the  10  orbits  deter- 
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mined  by  Hiller  during  this  time  interval;  all  the  observations  he  used,  and  others,  are 
included  in  che  new  determinations.  These  18  orbits  are  supplemented  by  13  orbits  deter¬ 
mined  by  Hiller  (with  che  same  computer  program)  for  dates  before  April  3  and  after 
August  17.  Thus  the  analysis  of  the  14th-order  resonance  utilizes  31  orbits,  for  dates 
between  1971  February  21  and  September  19.  The  variations  in  both  inclination  and  eccen¬ 
tricity  are  analysed,  but  the  inclination  proves  to  be  the  more  satisfactory  of  the  two. 

2  THE  18  NEW  ORBITS 

Observations  of  1970-97B  were  available  from  a  wide  variety  of  sources:  993  obser¬ 
vations  were  used  in  determining  the  new  orbits,  58Z  being  from  the  US  Navy  Navspasur 
system,  33Z  from  volunteer  visual  observers,  3Z  from  the  Hewitt  cameras  at  Malvern  and 
Edinburgh,  5Z  from  the  theodolite  at  Jokioinen,  Finland,  and  1Z  from  the  kinetheodolite 
at  che  South  African  Astronomical  Observatory.  Of  these  observations,  those  from  the 
Hewitt  cameras  were  by  far  the  most  accurate,  but  the  last  two  groups  were  also  most 
valuable  in  providing  a  wide  spread  in  latitude  to  improve  the  reliability  of  the  orbits. 
Residuals  of  the  observations  were  very  similar  to  those  of  Hiller2  and  are  not  given 
again. 

The  18  new  sets  of  orbital  elements  are  listed  in  Table  1,  with  their  standard 
deviations.  The  raw  values  of  inclination  are  plotted  in  Fig  I.  The  main  perturbation 
in  inclination  for  a  high-drag  satellite  is  usually  a  decrease  due  to  che  effects  of 
atmospheric  rotation;  but  from  Pig  1  it  appears  that  the  resonance  probably  causes 
a  considerable  increase  in  inclination,  which  exceeds  the  decrease  due  to  atmospheric 
rotation. 


The  accuracy  of  Che  IS  orbits  may  be  assessed  from  Che  average  standard  deviations 
of  orbital  elements.  The  eccentricity  has  an  average  sd  of  0.000012,  which  is  equivalent 
to  a  distance  of  90  m  radially.  The  average  across-track  error  is  given  approximately 
by  the  average  of  the  standard  deviations  in  inclination  and  right  ascension,  which  is 
0.0013°,  equivalent  to  160  m.  For  an  orbit  of  such  high  drag  and  high  inclination  these 
are  excellent  accuracies. 

The  values  of  perigee  distance  a(!  -  e)  for  the  18  new  orbits  and  the  13  orbits 
of  Hiller  are  plotted  in  Fig  2.  The  standard  deviations  are  smaller  than  the  radius  of 
the  circles.  Ppima  facie,  the  values  look  reliable  and  the  curve  joining  them  is  free  of 
ugly  discontinuities  in  slope.  However,  there  is  one  peculiarity.  The  main  oscillation 
is  attributable  to  the  odd  zonal  harmonics  in  the  geopotential,  but  the  expected  ampli- 
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tude  of  this  oscillation  is  6.5  km,  whereas  the  decrease  from  the  first  peak  to  the  mini¬ 
mum  is  nearly  18  km,  instead  of  the  expected  13  km.  This  gives  forewarning  of  likely 
problems  in  interpreting  variations  in  eccentricity,  because  neither  resonance  nor  air 
drag  would  be  expected  to  alter  the  perigee  distance  by  as  much  as  5  km.  (The  shifting 
of  the  minimum  away  from  the  w  ■  90°  point  could  well  be  due  to  resonance  effects, 
however. ) 

3  ANALYSIS  OF  INCLINATION 

3. 1  Removal  of  perturbations 

To  obtain  values  of  lumped  14th-order  harmonics,  the  values  of  inclination  are 
fitted  with  the  aid  of  the  THROE  computer  program^,  after  removing  the  effects  of  various 
other  perturbations.  These  are: 

(a)  the  lunisolar  gravitational  perturbations; 

(b)  perturbations  due  to  zonal  harmonics  in  the  geopotential; 

(c)  the  J2  2 perturbation  incorporated  in  PROP; 

(d)  perturbations  due  to  atmospheric  rotation; 

(e)  the  effect  of  the  precession  of  the  Earth's  axis,  and 

(f)  the  effect  of  any  lumped  harmonics  of  order  14  not  included  in  the  THROE  model. 

Ocher  perturbations,  such  as  solar  radiation  pressure  and  tidal  effects,  are  treated  as 
negligible. 

The  perturbations  (a)  and  (b)  above  have  been  calculated  by  running  the  PROD  computer 
program  with  one-day  integration  intervals,  and  the  Jj  2  perturbations  are  available  with 
the  PROP  output.  These  perturbations  are  fairly  large, ’reaching  0.0023°  for  (a)  plus  (b), 
and  0.0017  for  (c).  The  values  of  inclination  from  the  31  orbits  after  removal  of  these 
three  perturbacions  are  shown  in  Fig  3. 

The  perturbations  (d)  and  (e)  are  removed  within  THROE.  The  existence  of  the  per¬ 
turbation  (f)  depends  on  the  model  adopted  for  the  THROE  computer  run,  which  is  discussed 
below. 

The  terms  to  be  included  in  the  resonance  model 

The  perturbations  of  inclination  caused  by  14th-order  harmonics  in  the  geopotential 
express  themselves  through  terms  of  the  form 
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sin 


(Y®  -  qw)  , 


(1) 


where  m  is  the  argument  of  perigee,  and  $  is  the  resonance  angle,  given  by 


♦  -  M  +  w  +  14(0  -  v)  , 


(2) 


M  being  the  mean  anomaly,  0  the  right  ascension  of  the  node  and  v  the  sidereal  angle. 

The  parameters  y  and  q  are  integers,  and  explicit  formulae  for  the  magnitude  of  the 

perturbations  due  to  the  most  important  pairs  of  values  of  (y,q)  have  been  given  pre- 
7  8 

viously  ’  and  will  not  be  repeated  here.  These  formulae  depend  on  the  lumped  harmonics, 

g 

and  the  magnitudes  of  these  s>ay  be  very  roughly  assessed  by  assuming  that  the  individual 

-  ~5  2 

coefficient  of  degree  l  and  order  m  is  of  magnitude  10  /l  .  The  relative  magnitudes 

of  the  main  (y,q)  terms  for  1970— 97B,  calculated  on  this  basis,  are  given  in  Table  2. 

Table  2 

Orders  of  magnitude  of  the  perturbations 
in  inclination  for  1970-97B, 
for  various  (Y.*qT 


Here  the  main  terms  are  likely  to  be  those  with  (y,q)  -  0,0),  (1,1)  and  (1,-1),  but 

because  of  the  high  eccentricity  of  1970-97B  (e  *0.1),  the  q  «  2  terms,  which  are  of 

2 

order  e  ,  are  unusually  large  and  cannot  be  ignored.  However,  we  cannot  expect  to  deter¬ 
mine  five  pairs  of  lumped  harmonics  with  only  31  orbits  to  analyse.  So  the  THROE  fitting 
is  made  with  the  three  main  terms  already  specified,  leaving  the  (Y,q)  -  0,2)  and  (1,-2) 
terms  to  be  estimated  separately  and  removed.  Fortunately  the  removal  of  these  terms  could 
be  achieved  with  good  accuracy,  because  (a)  their  effect  was  quite  small,  the  total  change 
in  inclination  never  being  more  than  0.0005°,  and  (b)  the  value*  for  the  lumped  harmonics 
to  go  with  these  terms  were  fairly  well  known  beforehand10,  and  were  then  refined  by  the 
use  of  the  most  recent  solution  for  1 4th-order  harmonics 1 1 . 

So  THROE  was  run  with  (y,q)  -  (1,0),  (1,1)  and  (1,-1)  terms.  The  value  of  the  atmos¬ 
pheric  rotation  rate  A  within  THROE,  used  to  remove  the  perturbation  (d),  was  set  at 
0.9  rev/day  in  accordance  with  the  results  of  Ref  12.  Here  the  relevant  values  of  the 
parameters  controlling  A  were  as  follows:  perigee  height  near  230  las;  local  time  averaged 
over  three  complete  cycles2;  and  quite  a  strong  seasonal  bias  towards  sinner.  So  the  value 
of  A  is  given  by  the  curve  labelled  'average  local  time,  sunser'  in  Fig  12  of  Ref  12; 
whence  A  -  0.9  .  The  density  scale  height  H  at  a  height  }H  above  perigee  also  has  to  be 


specified  and  was  taken  as  37  km  in  accordance  with  the  COSPAR  International  Reference 
1 3 

Atmosphere  1972  for  the  appropriate  level  of  solar  activity,  which  gives  an  exospheric 
temperature  of  900  K. 

Although  A  varies  between  day  and  night,  the  constant  value  taken  here  should  not 
cause  any  errors  greater  than  the  observational  errors  of  i  ,  because  there  are  three 
day-to-night  cycles.  The  total  decrease  in  i  due  to  atmospheric  rotation  is  0.0120°, 
that  is  an  average  of  0.0020°  for  half  a  day-to-night  cycle.  Thus  if  A  is  in  error  by 
20Z  over  a  half  cycle,  the  error  introduced  would  be  0.0004°,  to  be  cancelled  by  an 
opposite  error  in  the  next  half  cycle,  if  the  mean  value  has  been  correctly  chosen. 

3 . 3  The  fitting  of  the  variation  of  inclination  with  THROE 

The  31  values  of  inclination,  after  taking  account  of  perturbations  as  already 
specified,  were  fitted  with  the  aid  of  the  THROE  computer  program  with  (y,q)  ■  (1,0), 

0,1)  and  (1,-1).  The  variations  of  the  resonance  angle  ♦  ,  and  of  i  ,  are  shown  in 
Fig  4. 

The  fitting  was  quite  satisfactory.  In  the  initial  run  the  measure  of  fit  e  was 
1.57  and  after  relaxation  of  the  accuracies  of  10  points,  the  value  of  c  was  reduced 
to  0.77  and  the  following  values  for  the  six  lumped  harmonic  coefficients  were  obtained: 

9—®  *  ^  , 1 

10  C[4  -  -  32  ±  24  ,  10*S14  -  -  20  ±  33  ,  (3) 

9J»0  1,0 

10  C]4  -  -  67  ±  53  ,  I0*S14  -  12  t  50 

9_-1.2  -1,2 

10  CJ4  -  -  77  i  53  ,  10  S,4  -  -  67  ±  40 

These  lumped  harmonics  are  linear  sums  of  individual  14th-order  harmonic  coefficients; 
explicit  expressions  are  given  in  section  3.4. 

The  fitting  of  the  curve  to  the  points  is  shown  in  Fig  5.  In  most  respects  the 
fitting  is  good.  First,  there  is  e  clear  increase,  of  0.010°,  between  the  early  values 
and  those  after  resonance:  this  is  eight  times  the  average  standard  deviation  of  the 
values,  so  there  is  a  well-defined  change  in  inclination  at  resonance,  which  should  be 
well  determined.  Second,  the  first  six  values  are  consistent  with  a  constant  value  of 
inclination,  and  so  are  the  last  eleven;  this  suggests  that  the  perturbations  have  been 
removed  successfully,  and  in  particular  that  the  right  value  of  A  has  been  used.  Third, 
although  10  values  were  relaxed  in  accuracy,  most  of  the  relaxations  might  be  called 
'cosmetic',  only  three  of  the  values  (at  MJD  41085,  41093  and  41141)  being  more  than 
twice  their  original  standard  deviations  distant  from  the  curve.  When  the  character  of 
the  fitting  became  clear,  several  of  the  orbits  were  re-run,  to  see  whether  omission  of 
some  observations,  or  relaxation  of  their  accuracy,  would  lead  to  values  of  inclination 
which  were  of  smaller  standard  deviation  and  also  closer  to  the  curve*.  The  values  given 
in  Table  1  record  the  final  results  of  this  process.  It  was  found  that  the  point  at 


MJD  41085  strongly  influenced  the  results  when  used  with  its  original  standard  deviation, 
and  in  view  of  its  'nonconformity',  we  decided  to  nullify  its  effect  by  increasing  the 
standard  deviation  by  a  factor  of  10. 

Although  the  change  in  inclination  -  almost  exactly  0.01°  -  is  quite  accurately 

defined  by  the  fitted  curve,  our  attempt  to  evaluate  six  lumped  coefficients  from  such 

sparse  data  is  obviously  liable  to  be  over-ambitious  -  there  are  only  13  points  in  the 

region  where  the  main  resonant  variation  occurs.  So  it  is  not  surprising  that  several 

of  the  values  (3)  and  (4)  are  numerically  indeterminate.  Unfortunately  none  of  the 

pairs  of  values  is  negligible:  but  since  the  values  of  the  (y,q)  “  (1,0)  coefficients  are 

quite  well  know1®’11,  we  calculated  that  effect,  corrected  the  values  of  inclination 

accordingly,  and  then  made  another  THROE  fitting  with  (v,q)  "  (1,1)  and  (1,-1)  only. 

9-0 , I  9-0,1 

The  (r,q)  ■  (1,0)  terms  (with  10  CJ4  ■  -6.5  ,  and  10  ■  -  15.7)  produce  an  increase 

in  inclination  of  0.0021°,  which  is  20Z  of  the  total  increase  in  Fig  5,  and  so  the  fitting 
with  (y,q)  «  (1,1)  and  (1,-1),  shown  in  Fig  6,  has  an  increase  202  smaller  than  in  Fig  3. 
But  Fig  6  is  otherwise  similar  to  Fig  5,  thus  suggesting  that  the  fitting  is  quite  stable. 

The  values  for  the  lumped  harmonics  obtained  from  the  THROE  fitting  of  Fig  6  are: 


o_  ■  ,u 

10*C,.  -  -  49  ±  48  , 

1 4 

9_“ 1 » ^ 

10  C|4  -  -  97  ±  48  , 

with  c  »  0.76.  Although  the  standard  deviations  are  slightly  lower  in  equations  (5)  than 
in  equations  (4),  the  earlier  solution  seems  better  in  principle,  because  the  removal  of 
the  important  (y,q)  “  (1,0)  terms  is  rather  drastic  surgery.  The  significant  feature  is 
the  complete  consistency  between  the  sets  of  values  (4)  and  (5).  This  indicates  that  the 
fitting  of  six  coefficients,  although  a  dangerous  procedure,  was  safely  accomplished. 

We  therefore  opt  for  the  solution  given  in  equations  (3)  and  (4) . 


oJ.O 

10  S,.  -  22  *  44 

1 4 


9-~  * ’ 2 
10  S, , 

1  4 


70  t  29 


(5) 


3.4  The  equations  for  the  lumped  harmonics 

The  six  lumped  harmonic  coefficients  appearing  in  equations  (3)  and  (4)  are  linear 
sums  of  individual  14th-order  coefficients,  and  ^  ,  and  are  given  by  the 

following  equations,  in  which  the  numerical  values  are  obtained  with  the  aid  of  the  PROF 
computer  program: 

.0,1 

C14  "  C 1 5 , 1 4  *  °*472£17,I4  +  °-057519,l4  *  0’152521,14  “  °*,98£23,14  “  °',52£25,14  > 


(6) 


>■> 

“14 


C1 4, 1 4  ’  °-830C,6.l4  -  °*813C18,14  “  °'343£20,14  *  °-083£22,14  *  °*298£24,14 


*  0.308C.,  . .  *  0. 199C _  . , 


-  0.369C 


24,14  • 


'14 


'14,14 


+  0. 182C 


16,14 


0.209C 


18, 14 


20,14 


0.358C,,  )4  -  0.246C 


(8) 


The  equations  for  the  S-coef f icients  are  obtained  on  replacing  C  by  S  throughout. 

The  right-hand  sides  of  equations  (6)  to  (8)  are  terminated  when  the  expected  contri- 
bution  of  Ci  4  ,  which  is  assumed  to  be  of  order  10  /i  ,  falls  to  less  than  52  of  che 
main  term  (and  remains  below  52  for  higher-degree  terms). 

Equations  (6)  to  (8),  with  the  numerical  values  for  the  lumped  harmonics  from 
aquations  (3)  and  (4),  have  been  used  in  the  recent  solution  for  individual  harmonic 
coefficients  of  14th  order11.  For  five  of  the  six  lumped  harmonics  the  new  values  for 
the  individual  coefficients  fit  the  lumped  harmonics  to  within  1.1  standard  deviations. 
(The  exception  was  S^’  ,  where  the  residual  was  1.4  standard  deviations.)  We  may 
therefore  conclude  that  the  standard  deviations  in  equations  (3)  and  (4)  give  every 
indication  of  being  realistic. 

The  accuracy  in  geoid  height  implied  by  the  standard  deviations  in  equations  (3) 
and  (4)  may  be  roughly  assessed  by  multiplying  by  the  Earth  radius,  6378  km  and  dividing 
by  the  root  sum  of  squares  of  the  numerical  coefficients  in  each  of  equat'  ns  (6)  to  (8). 
On  this  basis  the  average  accuracy  in  geoid  height  implied  by  the  three  sets  of  lumped 

harmonics  in  equations  (3)  and  (4)  is  15  cm  for  (C,S)^*'  ,  20  cm  for  (C,S)J’^  and 

-  -  -1  2  4  4 
25  cm  for  (C,S)  ’  .  These  accuracies  are  not  nearly  so  good  as  for  the  slower  reson- 
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ances,  such  as  that  of  1970-47B  ,  but  are  still  useful  in  providing  constraints  for 
inclination  74°. 

4  ANALYSIS  OF  ECCENTRICITY 

The  eccentricity  decreases  from  0.095  at  MJD  41003  to  0.079  at  MJD  41213,  mainly  as 
a  result  of  the  effect  of  air  drag,  and  there  seems  little  hope  of  removing  this  huge 

perturbation,  equivalent  to  115  km,  with  sufficient  accuracy  to  analyse  the  resonant  vari¬ 

ation  in  eccentricity,  which  is  unlikely  to  exceed  1  km.  This  expectation  was  confirmed 
by  trial  runs:  when  the  raw  values  of  e  were  used  and  corrected  for  air  drag  within 

THROE,  several  attempts  all  led  to  values  of  e  near  30. 

So  a  different  approach  was  adopted.  The  perturbations  Ae  in  e  due  to  odd  zonal 
harmonics  and  lunisolar  perturbations  were  calculated  by  PROD  at  1-day  intervals,  and  che 
resulting  values  of  aAe  were  added  to  the  perigee  distance  a(l  —  e)  to  give  values  of 
perigee  distance  adjusted  for  the  effects  of  zonal  harmonics  and  lunisolar  perturbations. 
The  values  of  this  adjusted  perigee  distance,  Q  ,  are  plotted  in  Fig  7  and  should  show 
the  variations  due  to  air  drag  and  resonance. 

The  effect  of  air  drag  is  to  produce  a  slow  decrease  in  Q  ,  given  by  equation  (6) 
of  Ref  14,  which  may  here  be  written  with  sufficient  accuracy  as 

HM,  , .  , 

■  -  3K7f  {r^-f  *  °<0.04>}  , 


d£ 

dt 


(9) 
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where  H  is  the  density  scale  height.  For  1970-97B  the  average  value  of  Mj  between 

2 

MJD  41003  and  41213  is  0.65  deg/day  ,  so  that,  on  taking  e  =  0.088  and 
M  -5100  deg/day,  we  obtain  the  mean  dQ/dt  from  (9)  as 


5100  deg/day,  we  obtain  the  mean  dQ/dt  from  (9)  as 

=»  -  0.00040  H  da 


The  density  scale  height  has  to  be  evaluated  at  a  height  1 -5H  above  perigee,  and  from 
CZPA  1372  l3,  this  gives  H  =«  45  km  for  1970-97B  (with  exospheric  temperature  900  K) . 

Thus  Q  should  decrease  at  approximately  18  m/day,  giving  a  decrease  of  3.8  km  over  the 
210  days  of  the  analysis.  The  actual  decrease  is  about  5  Van,  which  is  not  very  different, 
but  the  form  of  the  variation  in  Fig  7  is  most  peculiar,  with  a  very  rapid  decrease  at 
about  60  m/day  for  the  first  60  days,  up  to  MJD  41060.  followed  by  a  much  slower  decrease 
at  about  14  m/day  up  to  MJD  41160,  and  then  an  even  slower  decrease  for  the  final  50  days. 

The  steep  slope  in  the  first  60  days,  if  caused  by  air  drag  alone,  would  indicate 
a  scale  height  of  about  150  km,  which  is  not  credible.  So  the  rapid  decrease  is  puzzling, 
and  the  possibility  of  propellants  escaping  from  the  rocket  during  its  first  months  in 
orbit  cannot  be  ruled  out.  The  first  six  points  in  Fig  7  have  therefore  been  ignored. 

Fortunately,  the  decrease  in  perigee  distance  over  the  main  resonance  period 
(MJD  41060-41160)  has  a  reasonable  slope,  and  we  removed  the  effects  of  air  drag  by 
assuming  that  drag  caused  a  linear  decrease  in  Q  ,  as  shown  by  the  unbroken  line  drawn 
in  Fig  7.  If  (4 r  )  denotes  the  decrease  in  perigee  distance,  relative  to  the  value 

P  at 

at  MJD  41061,  as  given  by  this  line,  we  may  define  a  revised  perigee  distance 
(a(l  -  e)] rev  given  by 

U(1  *  e),rev  ’  Q  +  (Arp)at  *  (11) 


and  this  quantity  would  reflect  the  variation  in  perigee  distance  produced  by  resonance 
and  smaller  perturbations,  if  the  atmospheric  perturbation  assumed  was  correct.  We 
then  obtain  revised  values  of  e  ,  namely 

e  »  1  -  —  [  a  ( 1  -  e)  ]  ,  (12) 

rev  j  '  rev  ’  v  ' 

where  a  is  taken  as  7235  km.  These  values  of  e  were  then  fitted  by  THROE,  with 
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drag  suppressed,  and  with  (y,q)  ■  (1,1)  and  (1,-1\  which  are  by  far  the  most  important 
(y,q)  terms  producing  variation  in  e  . 

The  first  THROE  run  fitted  quite  well,  apart  from  the  last  five  points;  so  these 
were  dropped.  Of  the  remaining  20  values,  two  had  large  residuals  and  were  relaxed,  one 
by  a  factor  of  2  and  the  other  by  a  factor  of  4.  The  following  values  of  the  lumped 
harmonic  coefficients  were  then  obtained,  the  value  of  e  being  1.40: 


-  49  ±  53 
19  ±  29 


,J.O 

,C14 

-  -  52 

i  40  , 

q.1,0 

10  S14 

9_->.2 

c1A 

-  -  53 

±  40  . 

9-” 1  * 2 

10  S,( 

The  fitting  of  the  curve  to  the  points  is  shown  in  Fig  8.  When  the  calculation  was 
repeated  with  the  decrease  due  to  air  drag  represented  by  the  broken  line  in  Fig  7, 
there  was  a  substantial  increase  in  e  ,  to  1.59,  and  in  the  standard  deviations.  A  solu¬ 
tion  with  lower  c  and  standard  deviations  would  no  doubt  be  obtained  by  using  a  line  of 
lower  slope  in  Fig  7;  but  this  scarcely  seems  justified,  since  the  unbroken  line  corre¬ 
sponds  to  a  scale  height  of  26  km,  which  is  already  at  the  lower  limit  of  credibility. 

3ecause  of  the  interaction  between  air  drag  and  the  resonance,  and  the  difficulty 
of  making  a  logical  choice  for  the  slope  of  the  line  in  Fig  7,  we  feel  that  the  values  (13) 
are  not  so  reliable  as  the  values  (4) ,  even  though  three  of  them  have  slightly  lower  stan¬ 
dard  deviations  than  the  values  (4).  However,  it  will  be  noticed  that  the  C-values  in  (4) 
and  (13),  are  consistent  to  within  half  the  sum  of  their  standard  deviations,  while  the 
S-values  only  differ  by  0.6  and  1.2  times  the  sum  of  their  standard  deviations.  This 
suggests  that  the  values  (13)  may  be  more  reliable  than  expected.  Also,  three  of  the  four 
values  (13)  are  within  one  standard  deviation  of  the  values  of  the  four  lumped  coefficients 

given  by  the  recent  solution  for  individual  harmonics''  -  namely  -3!  and  +31  for 
9--10  9--  -12 

10  (C,S).^  and  -40  and  -9  for  10  (C,S))4’  .  This  again  suggests  that  the  values  (13) 

are  reliable,  despite  the  rather  devious  methods  required  for  their  derivation. 

Although  the  values  (4)  have  already  been  used  in  the  recent  determination''  of 
individual  14th-order  harmonics,  because  the  results  (13)  were  not  then  available,  it 
would  probably  be  more  appropriate  in  any  future  determination  to  adopt  a  weighted  mean 
of  (4)  and  (13),  with  the  smaller  of  the  two  standard  deviations.  This  'compromise 
solution'  is 

9J.O 

10  C]4  -  -  58  ±  40  ,  10  sj4  "  -  18  ±  50 

9—  - 1 , 2  9— ~  *  > ^ 

10  Cu’  -  -  63  +  40  ,  10  sj4  “  -17  +  29 

These  values  are  all  within  1  standard  deviation  of  the  values  given  by  the  recent 
solution  for  individual  coefficients  '.  Their  average  standard  deviations  correspond 
to  accuracies  in  geoid  height  of  1 7  cm  for  (C,S)J’°  and  18  cm  for  (C,S)~4’2  . 

5  CONCLUSIONS 

The  orbit  of  I970-97B  has  been  determined  from  observations  for  31  epochs  during 
210  days  in  1971  when  the  satellite  was  close  to  14th-order  resonance.  These  31  orbits 
comprise  13  previously  determined  by  Hiller  and  18  new  ones.  Some  Hewitt  camera  and 
kinetheodolite  observations  were  available  and  the  orbital  accuracy  achieved  -  about 
90  m  radial  and  160  m  cross-track  -  was  good,  in  view  of  the  very  low  perigee  height 
(230  km). 

Because  of  the  high  drag,  the  orbit  passed  through  14th-order  resonance  rapidly, 
and  appreciable  resonance  effects  lasted  only  for  about  100  days,  covering  13  of  the 
31  orbits.  As  the  orbit  is  quite  eccentric  (e  =*  0.09),  six  lumped  harmonic  coefficients 
prove  to  be  of  importance  and  have  to  be  evaluated.  With  so  many  coefficients  in  a 
fitting  of  (effectively)  13  points,  the  values  obtained  are  bound  to  be  rather  inaccurate. 
However,  consistent  values  were  derived  from  analysis  of  both  inclination  and  eccentricity. 


and  there  is  every  sign  that  these  values,  as  given  in  equations  (3)  af.d  (14),  are 
reliable,  although  far  less  accurate  than  those  obtained  from  slower  resonances.  The 
orbital  inclination  of  1970-97B  was  74°,  and  the  lumped  harmonics  provide  useful  con¬ 
straints  for  this  inclination,  near  which  no  other  such  accurate  results  are  available. 
The  values  have  been  successfully  used  in  the  recent  determination  of  individual 
coefficients  of  order  14. 
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